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We investigate the apphcation of the singular value decomposition to compact-binary, 
gravitational-wave data-analysis. We find that the truncated singular value decomposition reduces 
the number of filters required to analyze a given region of parameter space of compact binary coales- 
cence waveforms by an order of magnitude with high reconstruction accuracy We also compute an 
analytic expression for the expected signal-loss due to the singular value decomposition truncation. 
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I. INTRODUCTION 

The coalescence of compact binaries composed of neu- 
tron stars and or black holes is a promising source of 
gravitational radiation for ground-based gravitational- 
wave (GW) detectors. The mass parameters of the GW 
signal are not known a priori. In order to detect GW 
from compact binary coalescence (CBC) events, a large 
number of filter templates are required to to probe the 
continuous component mass parameter space, (7711,7712), 
of possible CBC signals in the detector data to high 
fidelity [', "^J. Template waveforms are distributed in 
the space such that there is a small maximum loss of 
signal-to-noise ratio (SNR) (called the "minimal match") 
due to the mismatch between an arbitrary point in the 
mass parameter space and the nearest discrete point of 
the template bank. A standard choice for the minimal 
match is 97%, which, for a hexagonally tiled, flat, two- 
dimensional manifold, corresponds to neighboring tem- 
plates that have greater than 95% overlap. 

This redundancy implies that correlated calculations 
are required to filter the data with these templates. The 
singular value decomposition (SVD) can be used to elim- 
inate these correlations by producing orthogonal basis 
vectors that can be used for filtering and reconstructing 
the original template bank. 

This work will describe how to reduce the computa- 
tional redundancy in filtering the CBC signal parameter- 
space in order to more efficiently infer whether or not a 
GW is present. Specifically, we will explore a purely nu- 
merical technique using the SVD to reduce the number of 
templates required to search the data. We note that oth- 
ers have applied the use of SVD to GW data-analysis to 
analyze optimal GW burst detection [ , ] and coherent 
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networks of detectors [ ]. We also note that significant 
work has been done to analytically reduce the compu- 
tational filtering burden using interpolation for certain 
template waveforms [ , ]. 

This paper is organized as follows. First, we describe 
the framework for CBC filtering in the context of vec- 
tor inner products. Next we introduce the SVD as a 
way to reduce the number of filters required to approx- 
imately compute those inner products. We then derive 
an expression for the expected SNR loss in terms of the 
singular values. Finally, we demonstrate the application 
of this method to a set of CBC waveforms corresponding 
to binary neutron star (BNS) coalescences. 



II. METHOD 

A. Matched filtering 

CBC searches employ matched filtering as the first 
step in locating a GW signal [o]. The optimal filter- 
ing strategy weights both the detector output and tem- 
plate waveform by the inverse of the amplitude spectral 
density of the detector noise, a process called "whiten- 
ing". Representing both the whitened data and the a^^ 
whitened template waveform as discretely sampled time 
series, s = {si} and ha = {hai}^ respectively, the output 
of the matched filter at a specific point in time is given 
by the vector inner product 



Pa = hl- s. 



(1) 



In searches for GWs from CBC sources, the signals be- 
ing sought are chirping sinusoids with an unknown phase. 
The search over phase is accomplished through the use 
of complex-valued templates where ^h^ contains the 
cosine-like phase and ^h^ contains the sine-like phase. 
The filter output can be maximized over template phase 
by evaluating \pa\' 

In the absence of a GW signal, the whitened detector 
data consists only of noise, n = {n^}, and is a stationary, 



zero-mean, unit- variance, Gaussian random process, so 
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where ( ) denotes the ensemble average. When the tem- 
plate waveforms are normalized such that ^h^ ■ ^h^ = 
^ha • '^ha = 1, (2) yields 

(/J; • n) = , (3a) 

{{ma • nf) = {{^ha • nf) = 1 . (3b) 

When (3) is true, pa is called the SNR and indicates how 
likely it is that a signal is present in the data at that 
point in time [ ]. 

As explained in Sec. l^ h^ - h^' > 0.95 for adjacent 
templates. For those templates, pa and pa' differ by at 
most 5%. This suggests the existence of an approxima- 
tion scheme that would allow the SNRs to be computed 
to reasonable accuracy without explicitly evaluating all 
the template inner products. Next, we will look at how 
the truncated SVD can be used to replace the template 
bank with an approximate, lower-rank, orthogonal basis 
from which the SNRs can be reconstructed. 



B. Reducing the number of filters Avith truncated 
singular value decomposition 

The waveforms are parameterized by their component 
masses and we denote the a^^ template waveform of the 
M templates required to search a given parameter space 
as ha = {h{mi^ ^^2, ^i)}- Rather than filter the data with 
N = 2M real- valued filters {M complex- valued filters), 
we linearly combine the output of a basis set of fewer, 
real- valued, filters, iT^, to reproduce pa to the desired 
accuracy, p'^. The goal is to have 
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where A is the complex- valued reconstruction matrix we 
wish to find and the number of inner products is reduced 
from N to N' . In order to find the basis vectors, u^^ we 
use the SVD of the real- valued template matrix, H 



= {3?/ii,^/ii,3?/i2,^^2, ...,3?/^m,^^m} , (5) 

where p identifies rows of H and indexes the filter num- 
ber, and j identifies the columns of H and indexes sample 



points. In this definition, the row vectors J^( 



(2a-l) 



and 



H(^2a) are, respectively, the real and imaginary parts of 
the a^^ complex waveform, ha- An illustrative template 
matrix can be seen in Fig. 1. 
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FIG. 1. An example template matrix, H. Top: An illus- 
tration of how the input template time series is packed into 
the template matrix. Bottom: The matrix of the template 
time series where the y-axis indicates the template waveform 
and the x-axis indicates the time samples. It should be noted 
that these waveforms have been shortened and have not been 
whitened for illustrative purposes. 



The SVD factors a matrix such that [ , Sec. 14.4] 



N 



H^j = ^ v^j,aj,i 



(6) 



where v is an orthonormal matrix of reconstruction co- 
efficients whose columns, V/j^u, satisfy 



^VfiuVfiX = ^i.A , 



(7) 



a is a vector of singular values ranked in order of im- 
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FIG. 2. An example basis matrix, u. Top: An illustration of 
the resulting orthonormal basis vectors ordered from most to 
least important (bottom to top) in reconstructing H. Bottom: 
The matrix of basis waveforms produced by the SVD. The y- 
axis indexes the basis waveforms and the x-axis indicates time 
samples. It should be noted that these basis vectors have been 
computed from shortened, non- whitened template waveforms 
as mentioned in Fig. 1 purely for illustrative purposes. 



portance in reconstructing the H, and u is a matrix of 
orthonormal bases (e.g. an illustration can be found in 
Fig. 2) whose rows are basis vectors, u^^ satisfying 



H- 



fXjUjyj — d^j. 



(8) 



However, since a search for CBC signals only needs 
waveform accuracies of a few percent to be successful, it 



is possible to make an approximate reconstruction of H 



^fj^j ^ ^Ij •= X] ^/ii^cr^^^j , 



(9) 



u=l 



where N' < N. This reduces the number of rows of u 
used in the reconstruction. We create a new basis matrix 
u = {uiyj} = {ui^U2^ ...un'}^ where u indexes the filter 
number, j indexes sample points, and we have discarded 
the basis vectors that look least like the template wave- 
forms (i.e. with the lowest singular values). We can write 
(4) as 

Pa = (^(2a-l) - ^^(2a) j ' ^ 

= Yl (^(2a-l)r.Cr^ - ^V{2a)u(^u) (u^ ' s) , (10) 

where we have made use of the packing of H (5) and (9). 



C. Reconstruction accuracy 

As we are not reconstructing the original template 
waveforms exactly, there will be some inherent mismatch 
between H'^ and H^. We want to know the expected 
fractional SNR we will lose because of this difference. 

As stated previously, the inner product of a (normal- 
ized) template waveform, H^^ with itself is 



H^ • H^ 
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(11) 



where, in the second line, we have made use of the or- 
thogonality of basis vectors (8). A similar relation can 
be found for the inner product of the reconstructed wave- 



form, iJ^, with itself 



N' N 

#;■#; = ^«>2 = i_ ^ „2^^2^ (12) 

Because of the orthogonality of the basis vectors (8), the 
inner product between a template waveform, J?^, with a 
reconstructed waveform, H^^ is 



H^ - Hi - H'^ - Hj, - H'^ - Hi . 



(13) 



In addition, the two phases of the templates, which are 
packed adjacently in H (5), are orthogonal 



AT 



^(2/2-1) • ^(2/i) = ^ ^(2M-l)i.^(2M)iyCr^ = . (14) 



i/=l 



This implies that the inner product of the two phases of 
the approximate waveforms are given as 



N' 



^(2/i-l) ' ^(2/i) = X]^(2/i-l)iv^(2/i)z.cr^ 



u=l 



N 



Yl ^{2^-l)v^{2^)v(^l ' (15) 

v^N' + l 



The average fractional SNR loss, ^Pa.1 Pa.^ between a 
template waveform and the two phases of the same re- 
constructed waveform is given by 



bpa_ 

Pa 



:= 1 



\Pa 



(16) 



The following derives the mismatch in terms the of 
components we truncate from the SVD. First we com- 
pute these terms for a given signal waveform, s = 

3? (Ae^^/ic.), with phase, (j). The SNR from the exact 

waveform, |pq,(^)|, is given as 



wm 



^Wio^-A^i&'^Ky 



^ha ■ A^ (e'^f^ha) ^ 



1/2 



A, (17) 



in which we have used (11) and (14). The SNR from the 
approximate waveform, |p^|, is given as 



Km 



mr ■ A^ 



{e''f>K) ' 



9/J'„ • A^ U'^K 
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(18) 



We can expand (18) using the packing of H (5), (12), 



(13), and (15) to 

\Pam = A 



N 
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(19) 



Let us look at the higher order sums in (19). The sums 

^u=N'-\-i^fiu^u^ which are also found in (12), represent 

the power of vector H^ lost through the truncation of the 

SVD. These sums must be less than 1, Xl^^^v'+i ^fiu^l < 
1. However, since the objective is for the approximation 



to be such that 



Hn — H' \\ ~ 1%, we expect 



AT 



Yl ^M^^^ ^ ^ ' 

u=N'-\-l 



(20) 



and we can therefore drop terms that are higher than 
first order in these sums. Additionally, 



N 



Y ^f^^^f^'^^l 



u=N'-\-l 



< 



\ 



E ^>i)( E «;>^s)«i- (21) 



This means (19) is approximately 

N 
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(22) 



As physical signals will arrive in the detectors with 
random phases, we now average over the phase, ^, using 



1 r^"" 



(23) 



resulting in 



A, 



\Pa\ 



1 ^ 



/x=Ar'+i 



Substituting (24) in (16), we find the average fractional 
SNR loss for the a^^ template 
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1 ^ 

4 zZ (4«-i)m + 4«)mJ 
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The expected fractional SNR loss can be computed by 
averaging over the waveforms in the template bank using 



(24a) 


iU" 


: 1 1 1 : 




10-1 


^^^^^"\ 




10-2 


^\ ''■ 


(24b) 


^10-3 


\ \ \ 


tional 


"^ 10-4 


\ ■ 




10-5 


\ '-■ 


(25) 


10-6 


\ \ \ 




ir,-7 


1 1 \l 



M 

-T 



5pa 
Pa 



(26) 



Combining (25) with (26), remembering M = A/'/2, and 
using the orthogonality of reconstruction coefficients (7), 
we get 



1 ^ 

= — y 



27V ^ M- 

/i^AT' + l 



(27) 



It is not surprising that the expected fractional SNR 
loss is proportional to the square of the Frobenius norm 
of the truncation error of H 



IH-H'll 



E (^M. 
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E ^l- (28) 

v^N' + l 



The expected fractional SNR loss, (Sp/p), can be used 
as a threshold for deciding how many basis vectors to 
keep in the truncated S VD reconstruction of the template 
matrix. For detection purposes, we want (Sp/p) to be less 
than the minimal match of the template bank. 



III. APPLICATION TO COMPACT BINARY 

COALESCENCE GRAVITATIONAL- WAVE 

SIGNALS 

We apply the above procedure to BNS waveforms with 
chirp masses 1.125M0 < Mc < 1.240Mq and component 
masses IMq < mi^m2 < SM©. The number of tem- 
plates required to hexagonally cover this range in param- 
eters using a minimal match of 96.8% is M = 456, which 
implies a total number of filters N = 912. These non- 
spinning waveforms were produced to 3.5PN order [ ], 
sampled at 2048 Hz, up to the Nyquist frequency of 1024 
Hz. The last 10 seconds of each waveform, whitened with 
the initial LIGO amplitude spectral density, were used to 
construct H. 
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FIG. 3. The expected fractional SNR loss, {6p/p), given by 
(27) as a function of the number of basis vectors we retain (out 
of iV = 912). The region {Sp/p) > 10% should be ignored as 
the Taylor expansion of the fractional SNR loss in (22) is not 
valid in that regime. 



In Fig. 3, we plot (Sp/p) as a function of the number 
of basis vectors kept. If we require that {Sp/p) = 10~^, 
we find we can reduce the number of filters in the above 
template bank from N = 912 to N' = 118, about an 
order of magnitude reduction in the number of filters. 

In Fig. 4 we show how {Sp/p) compares to the actual 
distribution of Spa/pa^ where we have chosen random 
values of (j) for each template. We find it is a good mea- 
sure of the expected fractional loss of SNR. 

We have investigated how generic this reduction of fil- 
ters is for other regions of CBC mass parameter space 
(e.g., regions of parameter space with larger component 
masses), and find the reduction to be similar. We tested 
this by generating a template bank with a 96.8% minimal 
match, component masses between IMq and 34M0, and 
total mass below 35M0. We then ordered the templates 
by chirp mass, split the template bank up into patches 
of M = 456 templates, and computed the SVD for these 
patches. 

We can include larger portions of parameter space in 
the SVD by including more templates such that the num- 
ber of templates is smaller than the number of time- 
samples per template. However the compuational cost 
of the SVD of an A/" X L matrix with N < L grows as 
0{LN'^)^ thus including more templates nonlinear ly in- 
creases the cost. Another complication is that waveforms 
further apart in parameter space have smaller overlap. 
This will result in more basis vectors being required to 
reconstruct the waveforms to the same accuracy. There- 
fore, including larger portions of parameter space in a 
single SVD computation will result in diminishing re- 
turns for the computational cost. We propose to address 



IV. CONCLUSION 
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We have investigated how the SVD can be used to re- 
duce the number of filters needed when analyzing GW 
data for CBC signals. We have found the number of 
filters required to matched filter these template banks 
can be reduced by about an order of magnitude through 
truncating the SVD of these waveforms. This result dif- 
fers from other work that models CBC GW signals in 
approximate ways [ 1:^-14] by starting with an exact rep- 
resentation of the desired template family and producing 
a rigorous approximation with a tunable accuracy. 

We plan to explore several topics in future works. 
Among these are the derivation of a composite detec- 
tion statistic using only the SVD coefficients in order 
to minimize the computational costs associated with re- 
construction and the interpolation of signals not in the 
original template set. 



FIG. 4. Histogram of measured fractional SNR loss, 5pa/ pa^ 
where we have chosen a random values of cj) for each template. 
The mean value predicted by (27), shown as the dashed-black 
line, matches the measured mean shown in solid-red. 



this issue, as above, by breaking up the parameter space 
into patches for which we can independently compute the 
SVD, although how best to do this is beyond the scope 
of the present work. 
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